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We present the theory of damping of low-energy excitations of a trapped 
ON . Bose condensate at finite temperatures, where the damping is provided by 



the interaction of these excitations with the thermal excitations. We empha- 



a 

*^ • size the key role of stochastization in the behavior of the thermal excitations 

for damping in non-spherical traps. The damping rates of the lowest excita- 



tions, following from our theory, are in fair agreement with the data of recent 
JILA and MIT experiments. The damping of quasiclassical excitations is de- 
termined by the condensate boundary region, and the result for the damping 
rate is drastically different from that in a spatially homogeneous gas. 



After the discovery of Bose-Einstein condensation (BEC) (T|-§[, one of the major direc- 
tions in the physics of ultra-cold gases has been the investigation of collective many-body 
effects. Especially interesting is the behavior of low-energy collective excitations of a trapped 
condensate. The JILA |||| and MIT [||,[7[ experimental studies of the excitations related to 
shape oscillations of the condensate show that these excitations are damped and provide us 
with interesting results on the temperature dependence of the damping rates and frequency 
shifts. 

In this letter we develop the theory of damping of excitations of a trapped condensate 
in the Thomas-Fermi regime at finite temperatures, where the presence of a thermal com- 
ponent is important. We confine ourselves to the damping of low-energy excitations, i.e., 
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the excitations with energies E v <C fi, where \i is the chemical potential, and consider tem- 
peratures T Too (uj is the characteristic trap frequency) ranging almost up to the BEC 
transition temperature T c . Thus far, theoretical and numerical investigations of elementary 
excitations of trapped Bose-condensed gases predominantly remained on the mean field level 
|| |17| . The investigation of damping phenomena requires analysis beyond the ordinary mean 
field approach ||18|| . It should be emphasized that the damping of low-energy excitations in 



a trapped Bose-condensed gas differs fundamentally from the damping of Bogolyubov ex- 
citations in an infinitely large spatially homogeneous gas. In the latter case, characterized 
by a continuum of excitations, any given excitation can decay into two excitations of lower 
energy and momentum. This damping mechanism, first discussed by Beliaev for T = 
19fl and employed by Popov [^(J at finite temperatures, proves to be dominant at T <C //. 



In a trapped Bose-condensed gas the character of the discrete structure of the spectrum 
of low-energy excitations makes Beliaev damping impossible under conservation of energy. 
Therefore, irrespective of the relation between T and /i, the damping of excitations with 
energies 

E v 4Zn,T (1) 

has to be provided by their interaction with the thermal excitations. The damping mech- 
anism involves processes in which the low-energy excitation {u) and the thermal excitation 
(7) are annihilated (created) and another thermal excitation (7') is created (annihilated): 

i/ + 7<->y. (2) 

We will discuss the case where the thermal excitations 7, 7' are in the collisionless regime. 
Under the condition ([J) the energies E 1 of these excitations are much larger than the en- 
ergies E v of the low-energy excitations. Therefore, the damping mechanism governed by 
the processes fl2|) can be treated as Landau damping. For spatially homogeneous gases this 
mechanism was first discussed by Szepfalusy and Kondor 

It is worth noting that inside the condensate spatial region, at T < // the density of 
occupied states of thermal excitations peaks at the energies _E 7 ~T, whereas for T^>/i this 
happens at E^ ~ fi. As just the excitations with \i give the main contribution to the 
damping rate, the collective character of the thermal excitations remains important even at 

t > fi (cf. jnj). 



In a trapped Bose-condensed gas the damping of low-energy excitations is determined 
by the behavior of the wavefunctions and by the distribution of the level spacings of thermal 
excitations with energies _E 7 < fi, which depends on the trap symmetry. We emphasize that 
stochastization in the behavior of these thermal excitations plays a key role for damping in 
non-spherical traps. For quasiclassical [E v ^> Hu) low-energy excitations the main contri- 
bution to the damping rate Y v comes from the boundary region of the condensate, which 
makes Y v completely different from that in a spatially homogeneous gas. The damping of 
the lowest excitations [E v ~ UuS) is determined by the behaviour of the excitations in the 
entire condensate region. For this case the damping rates following from our theory are in 
fair agreement with the data of the JILA || and MIT |7j experiments. 

Elementary excitations of a Bose condensate trapped in an external potential V(r) are 
commonly defined within the Bogolyubov-De Gennes approach (see [24]]) on the basis of the 
grand canonical Hamiltonian 



2 



H 



~A+y(r)A(r)*W 
2m 2 



*(r) 



(3) 



assuming a point interaction between atoms, with U = ATcfi 2 a/m, m the atom mass, and a the 
(positive) scattering length. The field operator of atoms \& is represented as the sum of the 
above-condensate part iff' and the condensate wavefunction \&o = (^)- Omitting the terms 
proportional to \l/' 3 and and using the generalized Bogolyubov transformation ^'(r) = 
J2vbu u u{ r ) — blv*(r), where b u , bj, are annihilation and creation operators of elementary 
excitations, the Hamiltonian (H) is reduced to the diagonal form H — H +J2 U E u blb v , if the 



functions u v , v v satisfy the equations 

f-h 2 A 



2m 











~V U ~ 


)=E V 




Vv. 


Vv. 




Vv. 




-V u 



(4) 



The condensate wavefunction \l/o is determined by the well-known Gross-Pitaevskii equation. 
In the Thomas-Fermi regime, where \x ~ no m U (no m is the maximum condensate density) 
greatly exceeds hu, one has [25U2B|: ^ (r) = [( / u-^(r))/t/] 1 / 2 for/i>y(r) and zero otherwise. 
The low-energy excitations [E v <C fx) are localized in the condensate spatial region, and in a 
harmonic potential V(r) the functions f„ =u u ±v u can be written as p~6[] : 



2n (r)U 
E u 



±1/2/ v-1/2 



(5) 



where n (r) = |\I/o( r )| 2 is the condensate density, Zj = (2fi/mu 2 ) 1 ^ 2 the characteristic size of 
the condensate in the i-th direction, and u>i the i-ih trap frequency. 

Interaction between the excitations, caused by the terms proportional to \l/ /3 and \1/' 4 
Eq.(P|), leads to damping. Below we will assume the inequality 



m 



(T/n 0m U)(n 0m a 



3\l/3 



< 1 



(6) 



which is fulfilled up to T pa 0.9T C for the conditions of the JILA || and MIT |7j experiments. 
Just Eq.(|]) ensures that the contribution to the damping rate of low-energy excitations from 
the \1/' 3 terms is much larger than that from the \I/ /4 terms, and the damping is actually 
caused by the interaction of the low-energy excitations with thermal excitations through the 
condensate and governed by the processes (|2|). The interaction Hamiltonian responsible for 
these processes reads 



H iQt = U J d 3 r* *' t (* /t + 



(7) 



Second, under condition (|) the damping rate can be found within the first-order perturba- 
tion theory in H int : 



iNj^ \\ N _ N 
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J , h E~-Ey+E u +iO 



(8) 



where iV 7 = [exp (E 1 /T) — 1] 1 are equilibrium occupation numbers for the thermal excita- 
tions. The transition matrix element can be represented in the form 
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<y|Antk7>= | [sH^-iH^-H^'-Hy')} , (9) 

where H VT1 , = J <M (r)/^r)/-(r)./y* (r) and Hp = J d 3 r^o(r)f-(r)f+(r)f+*(r). 

Since energies of the thermal excitations Ey^>huj, these excitations are quasiclassical and, 
similarly to the spatially homogeneous case, one can write 



= (Ey/W^ + (no(r)U) 2 - n (r)U))^ f 7 (v) . 
Then, using Eqfl5|), from Eq.@ we obtain 

/ F f j \i/2 

< >={^j Jd 3 r<S> uy (r)f,(r)f;,(r), (10) 

where $^ 7 (r) = W^n/l^F^r), and 



A _ 2^ + K(r)^) 2 -n (r)^ + (no(r)[/)^ 
E 7A /R 2 + K(r)f/)2 



For the distribution of energy levels of the thermal excitations with a given set of quantum 
numbers 7 determined by the trap symmetry (in cylindrically symmetric traps 7 is the 
projection M of the orbital angular momentum on the symmetry axis) we will use the 
statistical Wigner-Dyson p7| , p8| approach which assumes ergodic behavior of the excitations. 
Then, the quantum spectrum of the thermal excitations is random and the sum in Eq.(H) can 
be replaced by the integral J dE-fdEy J2jV 9j9yRyj> in which g 7 (Ej) is the density of states 
for the excitations with a given set 7, and _R 7 y the level correlation function. In non-spherical 
harmonic traps g^E v 1 and, hence, R 7y > ~ 1. Then, putting N y — N y > = E u (dN 7 /dE 7 ) 
and writing (E 1 —Ey+E u +i0)~ 1 as the integral over time iJ^°dtexp{i(E^—Ey+E u +iO)t/h}, 
from Eqs. (Q), ( P~0| ) we obtain: 

r^^l^ReE/^^^/f^exp + i0) j / rfyr^(r)$;(r')J( 7 (r, r',t), (12) 

where the quantummechanical correlation function 

X 7 (r, r', t)^/^^ y e X p|i ^^'^ f^f^f^fy (v'). (13) 

In our calculation of the damping rate T v we will turn from the integration over the 
quantum states of the quasiclassical thermal excitations to the integration along the classical 
trajectories of motion of Bogolyubov-type quasiparticles in the trap. Following a general 
method (see |29|-|31~||), for if 7 (r, r't) we can write 
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where H(p, r)=y (p 2 /2m) 2 + 2no(r)U(p 2 /2m) is the Bogolyubov Hamiltonian, and r c i(i|r, p) 
the coordinate of the classical trajectory with initial momentum p and coordinate r. Then, 
Eq.fll^D is reduced to the form 

^wkH^J?^ {'f-)h^)Wt\r, P MK,-H(p, ,))g*. (15) 

We first consider temperatures T 3> /x, where the main contribution to the integral 
in Eq . fll5|) is provided by the thermal excitations with energies < /x. In this case the 
use of the statistical approach in non-spherical traps is justified by the fact that, as shown 
in the motion of corresponding classical Bogolyubov-type quasiparticles is strongly 
chaotic at energies of order /x. The characteristic values of p and t in Eq. (|T5| ) are of order 
(mno(r)f/) 1 / 2 and h/E u , respectevely. For quasiclassical low energy excitations (Ev^Tiuj) 
this time scale is much shorter than u~ l and important is only a small part of the classical 
trajectory, where the condensate density is practically constant and r c /(i|r, p) — r+yt, with 
v = dH/dp. Representing $^ 7 (r)$* 7 (r c z(i|r, p)) as |/ J/ (r)| 2 F 2 (r) cos (p v vt/fl), where p u is a 
classical momentum of the Bogolyubov-type quasiparticle, from Eq.(|T5D we obtain 

I\, = J d 3 r\f u (v)\ 2 T uh (v), (16) 

where T v h(r) is the damping rate of the excitation with energy E v in a spatially ho- 
mogeneous condensate of density n (r). For E u <C n (r)U we have []3~3f 



vh 



£? i/ 37T 3 ' 2 (rioa 3 ) 1 ' 2 T/ (AtiqUTi). A slight numericall difference from the Szefalusy-Kondor result 



21j| is due to the use of free particle Green functions in their calculation. 

For a trapped gas the result of integration in Eq.(^) drastically depends on the 
trapping geometry. For example, for cyllindrically symmetric harmonic traps |/;,(r)| 2 oc 
[E 2 + (no(r)[/) 2 ) -1//2 and strongly increases near the boundary of the condensate spatial 
region, where n U ^ E v . Just this region of distances determines the density of states gu 
and the damping rate T v . Here the contribution to T v from the thermal excitations with 
E 1 of order E v is important, and one should also include the Beliaev damping processes 
v <r-> 7 + 7'. Omitting the detailes of the calculation which rely on a generalized version of 
Eqs.(|l0|)-(|15D and will be published elsewhere, 

we present here the result for M = 0: 

r " =6 - 7 V7j n^wm- (17) 

For the lowest excitations (E v ~ hu) the characteristic values of p in Eq.([T5|) are of order 
(m/x) 1//2 , and the result of integration can be represented in the form 

IV = A v ^ T -{n, m a 3 fl\ (18) 

where A v is a numerical coefficient which depends on the form of the wavefunction of the low- 
energy excitation v. In contrast to the case of E v 3> hu, the calculation of A v requires a full 
knowledge of classical trajectories of (stochastic) motion of Bogolyubov-type quasiparticles 
in the spatially inhomogeneous Bose-condensed gas. 
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The criterion of the collisionless regime for the excitations with energies Ey ~ jj, assumes 
that their damping time T" 1 is much larger than the oscillation period in the trap cu -1 and, 
hence, the mean free path greatly exceeds the size of the condensate. From Eq.([H]) we find 
r^~(T/fo)(n 0m a 3 ) 1//2 and obtain the collisionless criterion 

(T/huj)(n 0m a 3 y/ 2 < 1. (19) 

Due to collective character of the excitations the criterion ([19]) is quite different from the 
Knudsen criterion in ordinary thermal samples. 

Remarkably, both Eq. ([[!]) and the assumption of stochastic behavior of thermal excita- 
tions with energies of order /i are well satisfied in the conditions of the JILA |5| and MIT 
@ experiments, where the temperature dependent damping of the lowest quadrupole ex- 
citations in cylindrically symmetric traps has been measured at temperatures significantly 
larger than fi. The JILA experiment 0, where the ratio of the axial to radial frequency 
(3 = uj z /ujp = v^8, concerns the damping of two quadrupole excitations: M = 2, E v = \f2uo pi 
and M = 0, E v = 1.8u p . Our numerical calculation of Eq. (|T5D , with W v from flBfl , gives 
A v m 7 for M = 2 and A v m 5 for M = 0. This leads to the damping rate r„(T) which is 
in agreement with the experimental data |J (see Fig.l). In the MIT experiment [|7|], where 
(3 = 0.08, the damping rate has been measured for the quadrupole excitation with M — 0, 
E v = 1.58u z . In this case we obtain A v w 10. The corresponding damping rate T V (T) 
monotonously increases with T and for the conditions of the MIT experiment |7[ ranges 
from 4 s _1 at T ~ 200 nK to 18 s _1 at T « 800 nK, which is in fair agreement with the 
preliminary experimental data. 

Importantly, under the condition (|i~9|) the damping rate T v of the low-energy excitations 
is much larger than the damping rate IV of the oscillations of the thermal cloud. This 
phenomenon was observed at JILA |J . One can easily find that for T ^> \x the damping 
rate Tt ^navr, where n~ (mT/2irh 2 ) 3 / 2 is the characteristic density of the thermal cloud, 
cr = 87ra 2 the elastic cross section, and vt ~ JT/m the thermal velocity. Accordingly, the 
ratio Tt/T u is just of order the lhs of Eq.(plJ). 

In spherically symmetric traps at any excitation energies one has a complete separation of 
variables, which means that the classical motion of Bogolyubov-type quasiparticles is regular. 
The excitations are characterized by the orbital angular momentum I and its projection M, 
and for given I, M the level spacing g~ x ~ Tnu can greatly exceed the interactions provided 
by the non-Bogolyubov Hamiltonian terms proportional to \l/' 3 and \l/ /4 . In such a situation 
the discrete structure of the energy spectrum of thermal excitations becomes important, and 
one can get nonlinear resonances instead of damping. On the other hand, stochastization of 
motion of thermal excitations can be provided by their interaction with each other or with 
the heat bath. In this case the damping rate r^flTB]) follows directly from Eq.(^) by using 
the Dyson relation for the level correlation function f^8j {g^E v ~ 1) and / 7 from the WKB 
analysis of Eq. (f|) . 

For T< fi the picture of damping of low-energy excitations changes, since T v will be deter- 
mined by the contribution of thermal excitations with energies E^T . In this case, the lower 
is the ratio T/fi the more questionable is the assumption of ergodic behavior of the thermal 
excitations. But, even if the stochastization is present, at T significantly lower than \i the 
temperature dependent damping of the lowest excitations will be rather small. For cylindri- 
cally symmetric traps from Eqs. (|l~2j), (|15|) one can find T v ~ (E v /h)(T / /x) 3 / 2 (n 0m a 3 ) 1 / 2 . 
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FIGURES 

FIG. 1. The damping rate T u versus T for the JILA trapping geometry. The solid (dashed) 
curve and boxes (triangles) correspond to our calculation and the experimental data [6] for the 
excitations with M = 2 (M = 0), respectevely. 
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